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Abstract 

In the present article we present exact solutions of the Dirac equation for elec- 
tric neutral particles with anomalous electric and magnetic moments. Using 
the algebraic method of separation of variables, the Dirac equation is sepa- 
rated in cartesian, cylindrical and spherical coordinates, and exact solutions 
are obtained in terms of special functions. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The experimental data on the existence of anomalous magnetic moment and electric 
dipole moment for the Dirac particle§ require for their complete understanding an exact 
description of the corresponding single-particle states, i.e., exact solutions of the Dirac 
equation with non-minimal interactions. In some physical situations the interaction between 
the anomalous magnetic moment and electric dipole moments with the electromagnetic 
tensor can be solved exactly. 

In fact, if for the electron and the muon (and apparently also for the r-lepton) the 
deviation of the magnetic moment from the Bohr's magneton and the value of the electric 
dipole moment are extremely small, on the other hand, for the neutrino {u^, t'^, i't) despite 
the magnetic moment is very small it is one of their nonzero characteristic parameters. Also, 
we have that for the hadrons the values of the normal and anomalous magnetic moments 
are of the same order of magnitude. For all the hadrons with large mean life it has been 
established, when experiments make it possible, the presence of electric dipole moment. 
For electrically neutral hadrons (neutron. A, H° hyperons) the interaction between the 
anomalous magnetic moment with the external field becomes significant. 

Among the above mentioned particles, the neutron plays a privileged position because of 
their physical applications. In fact, one of the most crucial problems of the present century 
is the security of the nuclear energy. One of the most important ways for analyzing the 
kinematic parameters of a nuclear reactor is the injection of neutrons and the subsequent 
measurement of the decay on time of the net population inside the reactor. The mathe- 
matical models for this problem are still primitive: the momentum is replaced by a time 
dependent delta function, and in general, the space distribution of neutrons in impulse is 
unknown for the experimenter. Since neutron beams are controlled by means of magnetic 
and electric fields, in order to understand the space distribution of neutrons it would be of 
help to know the corresponding wave equations, i.e, the solutions of the Dirac equation for 
the neutron. 
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In view of all the above mentioned, we will examine in the present article exact solutions 
of the Dirac equation for electrically neutral particles with non minimal interaction with 
an external electromagnetic field, that is, we analyze the interaction between anomalous 
magnetic and electric dipole moments with the electromagnetic tensor. 

The amount of articles devoted to the study of this problem is scarce, the reason is that, 
for solving the Dirac equation, as general rule, a complete separation of variables is required 
and, the inclusion of tensor field functions in the Dirac equation dramatically restricts the 
possibilities of separation of variablesi^ In this case, the possibilities are even more limited 
that those obtained for the Dirac equation minimally coupled to the electromagnetic fieldi, 
or in the presence of gravitational fieldsB This one is because the gravitational field goes 
into the Dirac equation in a geometric way, via the Lame's coefficients, and no additional 
matrices in the equation if we choose to work in a diagonal tetrad gauge. Regarding the 
vector electromagnetic fields, the inclusion of them in the Dirac equation via the generalized 
momentum (minimal coupling), does not conduce to the apparition of new matrices, and! 
the presence of tensor fields terms of the form g'y'^'y^'Tmn introduces in the Dirac equation 
a new functional dependence. The difficulties that arise are those we find when we include 
non geometrized fields like scalar, pseudo-scalar or pseudo-vector fields.SSi 

Now, we proceed to mention the most relevant exact solutions of the Dirac equation with 
anomalous moment reported in the literature, among them we have the constant magnetic 
field problem solved by Strocchiil and by Ternov et a0, a generalization of the Volkov's 
problem was found by Ternov et a and more recently this problem has been revisited 
by Barutlli. The problem in spherical coordinates for the central field has been studied 
by Ternov0 and BarutEUli, a good review on exact solutions of the Dirac equation with 
anomalous interaction can be found in0. 

In the present article we will discuss this problem from a unified point of view based on 
the algebraic method of separation of variablesi0, which consists in reducing the original 
Dirac equation to a sum of two commuting differential operators as follows 
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{Ki + K2}$ = 0, [KiK2]_ = 



(1.1) 



where each operator depends on a different set of variables. Then, we reduce the problem 



iterating this process, we are able to reduce a system coupled partial differential equations to 
four systems of ordinary differential equations. Regretfully, not always a complete separation 
is possible because of the structure of the metric functions and the form of the external fields. 

The article is structured as follows, in Sec II, the Dirac equation with anomalous electric 
and magnetic moments is separated and some exact solutions are exhibited. In Sec. Ill, 
the Dirac equation is separated in cylindrical and spherical coordinates, also, some exact 
solutions in terms of special functions are presented. In Sec. IV, we present a particular 
configuration where Dirac equation is separable and exact solutions are obtained. Finally, 
in the conclusions, we summarize and discuss our results. 



In this section we proceed to separate variables and to find exact solutions of the Dirac 
equation with anomalous electric and magnetic moments when this interaction can be writ- 
ten in Cartesian coordinates. Then, for a neutral particle interacting with a magnetic field, 
Dirac equation in Cartesian coordinates reads 



of solving the original equation to obtain solutions for Ki$ = A<I> = — where A is a 
constant of separation. Now, applying the scheme ( |1.1|) for Ki$ = A$ and K2$+A$ = and 



II. CARTESIAN COORDINATES 




(2.1) 



Applying the algebraic method of separation of 




we can write Eq. ( p.l|) as 



follows 




(2.2) 



with, 




(2.3) 
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(2.4) 



and 

where the magnetic vector can be written as, 

H, = H,{x,y) + H,{z,t) 



(2.5) 



(2.6) 



applying the condition V .H = 0, on (2^) we find that = Hk{t), and if we fix our attention 
to time independent fields this term can be omitted. In this way, A$ = —K^t^ becomes in an 
algebraic equation that establishes the relationship between the different spinor components 
of $ and and the value of the constant of separation A , which satisfies the relation 

= - E\ 

Choosing to work in the following representation for the gamma matrices 

/ 1 \ / \ / \ / a3 

f= ^ ,y= ^ ,7^= 3 1, (2.7) 

VO -1 / Via^ / VfT^ / Va^ 



we have that, 



[^^di-^^d2 + m^^^'' -gH,{x)-\]^ = 



takes the form 



{d^-gH^{x)-\)^i^2-{'m + iky)^i^^ = 

{d^ + gH,{x) + X)<^4;i-{m-iky)^i,2 = 
on the other hand, from A$ = —Kzt^ we find. 







( \ 


$2 




{tX - E)/k,<^, 


$3 




{tX - E)/k,^^ 






[ J 



(2.8) 



(2.9) 
(2.10) 



(2.11) 



the coupled system of equations ( |2.9| ) and ( p. 101) can be solved in terms of special functions 
for some values of the magnetic field Hz{x). Among them we have, 



(a) = (3, (6) = Px, (c) = (3/x, (d) = /?exp(r?x), 



(2.12) 



first, let us consider the simplest case of a constant magnetic field (a) Hz = p. Then, the 
solution of the system reads: 



$1 = ci exp(w'm2 + P + {g(3 + X)^x) + cs exp{-Jm^ + + {g(3 + X^x) (2.13) 



$2 = Ci 



-C2 



m2 + P + {gH -Xf-gH + X 



m + iky 



exp{Jm'^ + ky + {gH — A)^^) 



(2.14) 



m2 + A;2 + {gH - Xf + gH + X 
m + iky 



exp{Jm'^ + ^« + {9H — Xyx) 



{b) Hz = fix For this linear magnetic field we have that the system of equations (p.9|)(p.l0|) 
can be reduced to the following two second order differential equations 



[{d,, + gf3x + X){d, - g(3x - A) - {m' + k')]^^ = 



(2.15) 



[{d, - g(3x - X){d^ + g(3x + A) - {m' + A;^)]<l>4 = 



(2.16) 



after making the change of variables 



1/2 



gPx + X= \ — \ y 



we have that the equations ( ^.15] ) and ( |2.16| ) take the form 



^-}<l>i = 



2gP 



(2.17) 



(2.18) 



d^ y^ -gP + m^ + kl _ 



(2.19) 



the solution of the parabolic cylinder equation ( |2.18| ) can be written in terms of confluent 
hypergeometric functions M(a, 6, z )§ 
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$, = c,exp(-yV4)M(i + !^^, i ^) + c,,exp(V/4)M(l + !^^, | ^) 

(2.20) 

using (^^ ) and the recurrence relation 

(6 - l)M(a - 1, 6 - 1, 2) = (6 - 1 - z)M(a, 6, z) + zM\a, 6, z) (2.21) 



we obtain 



5$4 = Ciyexp(V/4)M(l + !^^, I ^) + c,exp(V/4)M(!^^, 1 ^) (2.22) 
where the constant 5 reads, 

5 = (2.23) 

The third case corresponds to the magnetic field (c) = (3/x. Substituting this expression 
into (|2]|) and (|2.10|) we obtain 



^ + M^-*M±i) -*;-A^-4*. = (2.24) 



the solutions of ( |2.24| ) and ( [^.25D also can be expressed in terms of confluent hypergeometric 



functions as follows, 
2 

(2.26) 



$4 = exp(-^) [a,y3^+^M{l + g(3-~k,2 + 2g(3, y) + y-^^Mi^gf] - ~k, -2g(3, y)} 



$i=exp(-|-){6o2/^^M(^/?-^, 2gl3, y) +Wy^~3^M{l-g(3-~k, 2-2gf3, y)} (2.27) 
where 

g(3X 



y = 2 JA2 + A;2 + m^x, k = , (2.28) 

^ ^ y^A2 + kl + m^ 
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and the coefficients ao, oi, bo and bi satisfy the following relations 

ai 2i 



1- 2(7/3) JA2 + P + m2 (2.29) 



bi ky — im V ^ 

Oq {iky — m) 



, — (2.30) 

^0 2{2gp+l)^\^ + k^ + m'^ 

finally, we have the exponential depending magnetic field given by (d) = P exp{rjx). Here 
we have that, after making the change of variables n = exp{7]x), the system of equations 
(PD ( ^ ) takes the form 



(iA;j^ + m)$4 - (rifid^ - (A + /?/i)$i = (2.31) 



{iky-m)<!>i + {r]fidf, +(A +/3/i)$4 = (2.32) 



the system of equations (p.31|) ( p. 321) is similar to the one obtained when we solve the Dirac 



equation for an electron in a magnetic field likeQ (d). Then we have that, 



$4 = exp( )ao/i " M(- H , h 1, ) 

rj T] rj rj rj 

+ exp( )ai/i " M( , hi, (2.33) 

^1 = mo , " . exp /i* . M - + ^ + 1, + 1, ^ 

ikj. — X , l3u. ,.feir,,,A ikx 2ikx 23a. 

- iai , . exp '< M ^ + 1, + 1, ^ 2.34 

fcy + 7712 7] 7] 7] 7] rj 



III. EXACT SOLUTIONS IN CYLINDRICAL AND SPHERICAL COORDINATES 

In this section we are going to find some exact solutions of the Dirac equation with 
anomalous moment in cylindrical and spherical coordinates. The Dirac equation in spheri- 
cal coordinates for a chargeless particle with anomalous electric moment, expressed in the 
diagonal tetrad gaugei'i, reads, 
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+ m + g-f^-r^E [> ^ = 



r sin 

substituting into equation ( p.lD the electric field given E by the expression E 
separating the angular variables from the radius and time, we have 



$ = 



smv 



7^7°$ = 



(3.1) 
e/r, and 

(3.2) 
(3.3) 



where is a constant of separation associated with the angular dependence. It is worth 



noting that the operator (3^) correspond to the one obtained in the separation of variables 
for the Dirac equation in the Schwarzschild metric0i3 when we work in the diagonal tetrad 
gauge. The spinor $ is related to by means of the expression, $ = 7^7'^\E'. 
Working in the standard representation of the gamma matrices. 
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i = 1,2,3 



(3.4) 



and substituting 



into (13.21) we obtain, 



d ^ i{ge + k) 
dr r 



dr 



i{ge + k) 



$1 - (m - E)$4 = 



+ im + E)$i = 



(3.5) 



(3.6) 



where the constant E is the eigenvalue of the operator idt .After substituting (|3.5|) into ( p.6| 
we obtain a second order differential equation given by the expression 

d"^ [ge + k)^ — i{ge + k 



dr' 



+ 



+ m-E' 



$1=0 



(3.7) 



whose solution0 can be expressed in terms of Zy{z) = aJy{z) + f3Ny{z) where Jv{z) and 
Nv{z) are the Bessel and Neumann functions respectively, and a and (3 are arbitrary con- 
stants. 



$1 = aor5Zi_.(^^+fc)(Vm2 - E^) 



(3.8) 



then, substituting ( p.8|) into (|3.5| ) we obtain, 

$4 



m 



-aor'2Z_ 



-i{ge+k) 



(3.9) 



here, it is worth noting that, given the matrix structure of eq. ( |3.2| ), the spinor components 
$1 and $4 are proportional to $2 and $3 respectively. 



Now we are going to solve Eq. (|3.1|) when a constant electric field E = e along the radial 
direction is present. In this case, the set of separated equations reads: 



d 



ik 



+ ige + 
ar r 



d . ik 
ige 

dr r 



$1 - (m - E)$4 = 



$4 + (m + E)<l>i = 



(3.10) 



(3.11) 



the solution of the coupled system of equations ( |3.10|) -( p.ll| ) can be written by the help of 
the Bessel functions, 



$1 = cr^ Zi_-^{\/m'^ — E'^r) 



m — E 



=$4 = cr2Z_i_j^(Vm2 — E'^r) 



(3.12) 



(3.13) 



The Dirac equation in cylindrical coordinates expressed in the diagonal tetrad gauge 
takes the form. 



(3.14) 



where we have included a term associated with the anomalous magnetic moment. Eq. (|3.14|) 
can be written as a sum of two commuting differential operators: 



(7SV«9. +7SV5o)$ = A<l>, ^ = 7V'^' 
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(3.15) 



(3.16) 



here, it is convenient to work in the following representation of the Dirac matrices, 
'a^ \ /a^ \ /a^ \ /O z 
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-aW V -a^ / \ -a^ / \ i 



then, after substituting ( |3.17| ) into ( p.l5| ) we arrive at 



r 



(3.18) 



-a'^dr + - + ima^ - gH + X)x = 
r 



where, taking into account ( |3.16|) we have 

e 



e 



V fc— a'-' / 



(3.19) 



(3.20) 



where = k1 — E'^. Then, substituting the explicit representation of the Pauli matrices 
into ( 3.181 ) we obtain that the equations governing the radial dependence of G are: 



{dr - - (m + i{\ - gH))Q2 = 

r 



(3.21) 



(dr + —)Q2 - (m - i(X - gH))Qi = 
r 



(3.22) 



Here, there are two simple cases for which the system ( 3.21| )-( p.22| ) can be solved exactly in 
terms of special functions. They are: a) the constant magnetic field H, and b) H = a/r. 
a) For this case we obtain the following second order differential equations 

/ + 1) 



dr'^ 



m'-{X- gHy $2 = 



(3.23) 



\dr^ r 
whose solutions take the formS 



^-^%^-m^-(A-^i/)2^ $1 = 



(3.24) 



$1 = cy/^Zk,.i/2{tJm^ + (A - gHyr 



(3.25) 
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$2 = c2r^'^Zk,+i/2{iJm? + (A - gHfr 



(3.26) 



where Za{z) is the general solution of the Bessel equation, and the constants ci and C2 are 
related as follows, 

m + i{\ — gH) 



Ci = IC2 



m2 + (A - gHf 



(3.27) 



b) Substituting H = a/r into (^)-(F^) we get 



{dr - y)ei - (m + zA - t^)Q2 = 



(3.28) 



{dr + ^)e2 - (m - zA + «^)ei = 



(3.29) 



the system of equations (p.28|) - (|3.29|) resembles the one obtained in solving the Dirac equa- 
tion with a Coulomb potential. Therefore it is possible to express the solution of (|3.28|) - (p.29|) 
by means of hypergeometric functions as follows: 



ei = i/l-^-e"i(Fi(p) + F2(p)) 
V m 



(3.30) 



ei = i/l + z-e"i(Fi(p)-F2(p)) 
V m 



where p = 2\/m? + A^r, and Fi and F2 read: 



F2{p) = Cop''M{-f 



ga\ 



+ A2 



, 27+1, p) 



(3.31) 



(3.32) 



Flip) = -co ^ P"M(7- 



(7aA 



+ 1, 27 + 1, p) 



(3.33) 



where Cq is a constant arbitrary . and 7 = y ^1 + S'^a^ 

Making use of the expression for the asymptotic behavior of the confluent hypergeometric 
functions for p —>■ oolli 



M{a,b,x) 



T{b - a 
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r a 



(3.34) 



we have that the solutions of ( p.2^ )-( p.29| ) are regular al infinity only if 

1 



or equivalently, 



(3.35) 



= kl - (3.36) 



where n =0, 1, 2, , and E is the energy. 

the expression ( p.36|) takes a very simple form when n = 0, 

E^ = kl + m'+[^-^^' (3.37) 

Now, we are going to solve the Dirac equation for a neutral particle with anomalous 
electric interaction. For an electric field along the z axis we have that the Dirac equation 
takes the form 

7°(9t + -i^dr + -^da + 7^(9^ + m + q^'^i^e] ^ = (3.38) 



Following the pairwise scheme of separation, we have that Eq. ( p.38|) can be reduced to a 
sum of two commuting first order differential operators 



[l^l^l^dr + ^ +gE - Aj $ = (3.39) 

(I'^do + 7°9, + m7=^7° + a) $ = (3.40) 

with $ = 7^7°\[' .Choosing to work in the Dirac matrices representation given by the ex- 
pression ( p.4|) we have that Eq. (|3.39| ) can be written as: 

(a^dr - yd^ - {gE - A) j 6 = (3.41) 



a' 



a^dr + {gE-\)\x = Q (3.42) 

r / 



with 
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e 



l(!B±^_3C) 



e 



(3.43) 



where the constant of separation A satisfies the relation: = /c^ + — E"^. From (|3.41|) 
we have that the components of G satisfy the relation: 



dr-^\ei = t{x-gE)e2 



(3.44) 



V + — ) 62 = -^(A - gE)Q, 
r / 



(3.45) 



the solution of this system reads 



ei = cr-2Zk,-i/2{i{X-gE)r) 



(3.46) 



82 = -cr2Zfc^+i/2(i(A - gE)r) 



(3.47) 



IV. A SOLVABLE EXAMPLE 

The following example could be regarded as surprising because of the big amount of 
functions appearing in the Dirac equation. This example seems to be contradictory since, 
following the remarks made at the introduction, the presence of additional matrices in the 
Dirac equation dramatically reduces the possibilities of separation of variables. But such 
a contradiction does not exist in the case to be presented. All the functions considered 
in the example, depend only on the radial variable r namely in the form 1/r, therefore 
their structure repeat the Lame function arising when we use spherical coordinates. This 
one explains the restriction in the election of the functions allowing complete separation of 
variables. 

Then, let us consider a Dirac particle in a Coulomb field with anomalous electric dipole 
interaction and also scalar and pseudoscalar interactions. At a first glance this problem 
seems artificial, because we do not relate the Coulomb term with the dipole one, but the 
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problem includes the possibility of a unified approach of a series of physical situations like 
the Hydrogen atom (Coulomb field), the non minimal electric dipole interaction, and con- 
finement model (scalar potential) and the quarkonium theoryil^ 

The Dirac equation expressed in the diagonal tetrad gauge takes the form 



1 



+ ±(y5, + ^9^)+7°(9*-^) + g7V- + C- + e7 



S 



^ = (4.1) 



here Ze'^/r is the Coulomb potential, qe/r is the electric dipole interaction, S is the scalar 
potential, P is the pseudoscalar potential, and ( and C, are constants. After separating the 
angular dependence, ( [4.1| ) takes the form 



(4.2) 



l^lWd^ + ^S^^^ - 7'7'e^(^, <^) 1 $ = fc$ 



(4.3) 



y sin-i? J 

where \1/ = 7^7°$ 

We will not solve the equation ( |4.ij| ) governing the angular dependence of the problem 
because this one is not possible without knowing the P('(9, ip) function. Obviously, the 
solution for the particular case P = will be expressed in terms of spherical harmonics like 
the Coulomb case. 



Choosing to work in the representation 



'ia^ 



cr 



7 



7 



iaV V 

equation ( |4.2| ) takes the form, 

k — eq 







-cr 



, 7 



a' 
-a^ 
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$1 = 



substituting ( [4.4| ) into ( |4.5| ) we arrive at 
\ d k — eq 



dr r 



61 +[-iE + m)- -{Ze' + C5) ) 02 = 



'0 1 
1 



(4.4) 



(4.5) 



(4.6) 



+ 62 + ({E - m) + -{Ze' -CS)]Q, = 

dr r \ \ r 



(4.7) 
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where 0i and 02 are the components of the spinor $i 

$1 



(4.8) 



Introducing the notation i( A; — eg) = M, m+E = A, m—E = B, (S+Ze"^ = a, Ze^ — ^S = /3, 
and rD = p we have that eq. ([4.6|) , (^4.7| ) take the form 



dp p ] ^ \ D p) ^ 



(4.9) 



d Ml ^ f B ^ 



(4.10) 



we shall look for solutions of the system ( [4 .91) , ( [4.10 ) in the form of power series: 

oo oo 

01 = E 02 = e-" P'^^'bv (4.11) 

v=0 v=0 

Substituting ( [4. Ill ) into ( [4.9| ), ( |4.1(J| ) and putting the coefficients of p*+"+i equal to zero we 
find 



[s + M)ao + pbo = 



(4.12) 



aao - (s - M)bo = 



(4.13) 



a-v^i + {s + v + M)a.u - j^^v-i + = 



bv-i - {s + V - M)by + j^(^v-i + abv = 



from (|4.12|) and ( |4.13 ) it follows that the parameter s takes the value 



(4.14) 



(4.15) 



(4.16) 



where we have considered a positive root in order to avoid a divergent solution at the origin. 
From ( 4.14 ) and ([4.15 ) we find the following relation between the coefficients a„ and b^ 
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{D{s + v + M)+ Ba} = {B{s + v - M) - Dp} 6„ (4.17) 



The series ( [4.11| ) will have a good behavior at infinity if they terminate for a finite value 



Putting ttN+i = bN+i = in ( |4.14| ) and ( [4.15| ) with oat 7^ and ^at 7^ we arrive at 



then, from (|4.17|) , (^4.18|) and the definitions of B and D we arrive at: 

EZe^ = {s + N)^w? - ^2 + mCS (4.19) 



The expression ( [4.19|) is the condition of quantization of the energy, and its structure is very 



similar to the one obtained for the hydrogen atom. From (|4.19|) can obtain some particular 
cases! a ) e = 0, S* = Hydrogen atom, b) e = 0, Z = Confinement of quarks by a scalar 
potential c) Z = 0, S* = 0. No minimal electric dipole interaction. 

Obviously, exact solutions of the Dirac equation containing such set of terms like ( |4.1| ) 
is not only of physical interest, but also is interesting from a mathematical point of view. 
Regretfully the amount of examples of this kind are relatively scarce and therefore not exact 
solutions are known from most of the problems. 



V. CONCLUDING REMARKS 

In this article we have obtained some exact solutions of the Dirac equation for electric 
neutral Dirac particle with anomalous magnetic moment. The results presented in this paper 
show the capabilities of the algebraic method of separation of variables provided a suitable 
representation for the gamma matrices. The separation of variables for other curvilinear 
system of coordinates was not possible because of the matrix character of the anomalous 
interaction in the Dirac equation. This peculiarity dramatically restricts the possibilities of 
finding exact solution in other geometric and physical configurations. Also, it would be of 
interest to find exact solutions when other interactions are present. 
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